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ABSTRACT 

Aims. Projected rotational velocities {vc sin i) have been estimated for 334 targets in the VLT-FLAMES Tarantula 
survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from 
approximately 09.5 to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, 
which is critical for understanding the evolution of massive stars. 

Methods. Projected rotational velocities were deduced from the Fourier transforms of spectral lines, with upper limits 
also being obtained from profile fitting. For the narrower lined stars, metal and non-diffuse helium lines were adopted, 
and for the broader lined stars, both non-diffuse and diffuse helium lines; the estimates obtained using the different sets 
of lines are in good agreement. The uncertainty in the mean estimates is typically 4% for most targets. The iterative 
deconvolution procedure of Lucy has been used to deduce the probability density distribution of the rotational velocities. 
Results. Projected rotational velocities range up to approximately 450 kms~^and show a bi-modal structure. This is 
also present in the inferred rotational velocity distribution with 25% of the sample having <We<100kms~^and the 
high velocity component having 250kms~^. There is no evidence from the spatial and radial velocity distributions 
of the two components that they represent either field and cluster populations or different episodes of star formation. 
Be-type stars have also been identified. 

Conclusions. The bi-modal rotational velocity distribution in our sample resembles that found for late-B and early-A 
type stars. While magnetic braking appears to be a possible mechanism for producing the low- velocity component, we 
can not rule out alternative explanations. 

Key words, stars: early-type - stars: massive - stars: rotation - stars: magnetic field - Magellanic Clouds and associa- 
tions: individual: Tarantula Nebula 



1. Introduction 

Rotation is thought to have a major influence on the evo- 
lution of massive stars. It may affect both how the star 
evolves dur ing hydrogen core burning (see, for example, 
iBrott et al.ll201ll: lPotter et al.ll2012bD an d the nature of any 
supernova explosion (s ee, for example IWooslev &: Heged 
[20061 lYoon et al.l [20061) . It can also provide insights into 
other phenomena such as star formation (see, for example. 
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IWolff et al."2007^ and how magnetic fields affect stellar evo- 
lution ((Potter ct al. 2012a). 

There have been several large studies of rotation in mas- 
sive main sequence stars both in our own Galaxy and in 
the Magellanic Clouds. For the former, one of the earlies t 
comprehensive samples was provided bv ISlettebakI ()1949[ ). 
As this sample was selected without any previous knowl- 
edge about the stellar spectra, it provides an unbiased sam- 
ple. Estimates of the projected rotational velocities of 123 
early-type stars were presented together with upper lim- 
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its for 63 stars with relatively narrow absorption lines. 
lAbt et all (|2002f ) investigated 1100 B-type stars selected 
from the bright star catalogue and this is one of the largest 
samples for which the biases are well understood. For late 
B-type stars they find evidence for a bi-modal distribu- 
tion but early-B type stars do not seem to show this bi- 
modalit y, although th e numb er of such stars in their sample 
is small. iHuang et al.l (|2010[ ) combined data from their ob- 
serving campaigns to provide projected rotational velocities 
for approximately 660 cluster stars and 470 field stars. They 
find that the average projected rotational velocity is signif- 
icantly lower for field stars, which they attribute to the fact 
that they are systematically more evolved and spun down 
than their cluster counterparts. 

For the Magellanic Clouds, there have been two investi- 
:ation s using ESO large programmes, viz. flVIartav an et al 
2006', '2007 1, and ref e rences there in) and IHunter et al 
2008b' ): .Evans et al.l (|2005l l2006l hereafter designated 
FLAMES-I). These observed predominantly field stars, 
while as discussed above for the Galactic studies, the be- 
haviour of cluster stars may be different. Also given the 
limited time cadence of the observations, it was not always 
possible to reliably distinguish between single and multiple 
systems. 

Here we present projected rotational velocity estimates 
for a sample of 334 early B- and late 0-typ e stars located in 
the 30 Doradus region (jEvans et al.ll201ll hereafter Paper 
I). The rotational velocities of the more massive 0-type 
stars will be discussed elsewhere (Ramirez-Agudelo et al. 
in preparation). Our sample has been selected as having no 
significant radial velocity variations, thereby reducing the 
contamination from multiple systems. We have also inferred 
the intrinsic rotational velocity distribution and discuss its 
bi-modal structure. 



2. Observations 

The spectroscopic data were obtained as part of the 
European Southern Observatory (ESO) large programme 
(|Evans et al.l 120111 ) using the Fibr e Large Array Multi- 
Element Spectrograph, FLAMES (jPasauini et al.l I2002D 
on the Very Large Telescope (VLT), primarily with the 
Giraffe spectrograph, but also supplemented with data from 
the Ultraviolet and Visual Echelle Spectrograph, UVES 
(jPekker et al.ll200"ol) . All observations have been based in 
and around the LMC open association 30 Doradus, with 
the core region R136 having a 15 arcsecond exclusion ra- 
dius as it is too densely populated for use of the Medusa 
fibres. This region has been observed with five pointings 
using the ARGUS integral field unit (IFU). The current 
work will use the observations obtained using the two bluer 
region Giraffe spectra to obtain estimates of the projected 
rotational velocity and the red spectra to search for Be- type 
emission. The observational configurations are summarized 
in Table m with details of the target selection, observations 
and data reduction being available in Paper L 

The sample is that defined in Paper I as potentially 
B-type stars and comprises approximately 540 targets. 
Ongoing spectral classification indicates that some targets 
are late-0 type stars with spectral types of 09.5 or 09.7. 
It is believed that our sample is effectively complete as re- 
gards the 09.7 and later spectral types but is incomplete 
for the 09.5 spectral type as approximately half these tar- 
gets were included only in the dataset of potentially 0-type 



stars. We return to this in Sect. 13.31 Additionally the mag- 
nitude cutoff employed in the survey implies that subject 
to spatial variations in extinction our sample should not 
include (near) zero-age main sequence stars with spectral 
ty pes later than a p proxim ately B3. 

iDunstall et al.l ()2012[ ) have undertaken a radial veloc- 
ity analysis for this dataset to identify binaries. They find 
37% of their sample show unambiguous velocity variations 
with an additional 6% having an uncertain binary status. 
We excluded the former from our analysis but included 
the latter, although we consider these results separately. 
Additionally thirty stars without significant radial velocity 
variations have been identified as supergiants and a prelim- 
in ary model atmosphe re analysis, similar to that employed 
by iFraser et al.l ()2010[ ) implied that they all had logarith- 
mic gravities between 2.5 and 3.2 dex. As these objects will 
normally have evolved from 0-type stars, they will be dis- 
cussed elsewhere (McEvoy et al, in preparation; Simon-Diaz 
et al., in preparation). 

Hence our sample of 334 stars (of which 45 are of un- 
certain status) is a subset of the Tarantula B-type sample 
that lie close to the main sequence. They show no convinc- 
ing evidence of any radial velocity variation typically at 
the 10 kms^^ level over six epochs spread over one year. 
Assuming a spectral type range of B3 to 09.5, leads to an 
effectiv e temperature range of approximately 20 000 to 34 

000 K (|Trundle et al.ll2007l) and a ma ss range of approxi- 
mately Vtol^'^I^^Brotreran [20111) for an LMC metal- 
licity; 

IDunstall et al.l ()2012[ ) have discussed completeness and 
estimate that up to 50% of the apparent B-type main se- 
quence sample could be undetected binaries. Indeed the 
failure to identify long period systems (typically P > 
1000 days) is an intrinsic limitation of the VLT-FLAMES 
Tarantula Survey (VFTS) given the one year baseline of 
the observations. Assuming that this analysis is correct, 
our sample could contain a significant fraction of binaries, 
albeit preferentially with long orbital period and/or low- 
mass companions. The light and spectrum from such bina- 
ries would be dominated by the main (B-type) component. 
Additionally the long orbital periods involved are unlikely 
to have led to any significant tidal interaction that would 
have modified the rotation of the objects. 

As there were no significant radial velocity variations in 
our sample, spectra could be combined without including 
any wavelength shifts. Several methods were adopted with 
the earlier analysis being based on co-added spectra with 
exposures containing cosmic rays or having a low signal-to- 
noise (SNR) being excluded. Later, all the exposures were 
normalised using low order polynomials. Normally an ex- 
tended wavelength region could be normalised with a single 
polynomial but for some exposures, there was significant 
variations in the continuum flux with wavelength (se Paper 

1 for details). In these cases individual lines or groups of 
lines were normalised. Spectra were then combined with a 
simple median filter or a weighted mean with cr-clipping; in 
both cases low SNR exposures were again excluded. Given 
the relatively small range in SNR for a given star and wave- 
length region, these two approaches gave essentially iden- 
tical spectra. Comparison of these data and of the earlier 
co-added spectra for five stars showed excellent agreement 
(at the 1-2% level) in the estimated mean projected ro- 
tational velocities indicating that the choice of reduction 
technique should not be a significant source of uncertainty. 
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Table 1. Summary of the spectroscopy used in the rota- 
tional velocity analysis. For some stars additional LR02 ob- 
servations were obtained when the initial observations were 
found to be of poor quality due to weather conditions. 



Instrument 


A-range 


R 




Exposures 




(A) 






(seconds) 


Giraffe LR02 


3960-4560 


7000 


6 X 


(2 X 1815) 


Giraffe LR03 


4505-5050 


8500 


3 X 


(2 X 1815) 


Giraffe HRI5N 


6445-6815 


16000 


2 X 


(2 X 2265) 



3. Projected rotational velocities 

3.1. Methodology 

Several methods are available for estimating projected ro- 
tational velocities (wesinz). Those most commonly used 
for OB-type stars include a direct measure of the full- 
width- half^mBxiimi^ spectral hues ( see, for ex- 
ampl e ISlettebak et~all Il975t lAljt et all l2002t IStrom et al.l 
20051) ■ profile fitting (see, for example, Rvans et al. 2002t 
Hunter et al.ll2008bD o r cross correlation (see, for example, 



Howart het al.lll997l) . ISimon-Diaz fc Herrerol jlQDl) have 
discussed the utility of the Fourier transform (FT) method 
for early type stars, as this should be able to separate the 
rotational broadening from other broadening mechanisms 
such as macroturbulence. Recently this approach has been 
widely used to study the different mechanisms contributing 
to the broade ning of spectral lines in early- type supergiants 
(iDufton et al.ll2006al: llefever et al.l 120071: iMarkova fc Puis! 



I2008t iSimon-Di'az et al.ll2010l: iFraser et al.ll2010[ )? 

We have adopted this Fourier transform method for 
all the stars in our sample. This approach relies on iden- 
tifying the first minimum in the Fourier transform for a 
spectral line, which is assumed to be the first zero in the 
Fourier transform of the rotational broadening profile with 
the other broadening mechanisms exhibiting either no min- 
ima or only minima at higher frequencies. Fig. [T] shows the 
He I line at 4026A and its Fourier Transform for two tar- 
gets, VFTS046 and VFTS699. Their Hen spectra are of 
similar strength (implying a spectral type near BO) and 
hence the difference in their line profile refiects their dif- 
ferent projected rotational velocities. Also shown are the 
rotational broadening profiles for a projected rotational ve- 
locity, corresponding to the first minimum in the Fourier 
Transform. As would be expected these agree well with the 
observations in the line centre but are narrower in the wings 
for VFTS046 due to the intrinsic Stark broadening of this 
diffuse line. 

As discussed bv ISimon-Diaz fc: Herre ro (2007), the first 
minimum in the Fourier Transform can be difficult to iden- 
tify in spectra with, for example, low SNRs or significant 
nebular contamination. Therefore we have also profile fit- 
ted (PF) all the spectral features, assuming that rotation 
dominates the line broadening i.e. we have used a simple 
rotational profile and assumed that the intrinsic and instru- 
mental profiles can be approximated by scaled (5-functions. 
The neglect of other broadening mechanisms will lead to 
these estimates being upper limits but they provide a use- 
ful constraint on the actual projected rotational velocities. 
We discuss these results further in 13.21 





4024 4026 4028 4030 4032 4034 4035 





4025 4030 4035 



Fig. 1. Profiles and Fourier transforms for the Hei line at 
4026A in VFTS046 and VFTS699. The theoretical pro- 
files correspond to the Vc sin i estimates from the Fourier 
Transform methodology. The position of the first mimium 
in the Fourier Transform is marked by a vertical line. For 
VFTS699, it has been necessary to remove some nebular 
emission near the line centre by linear interpolation over 
the affected wavelength region. 



Two different sets of spectral lines were adopted for 
analysis depending on the degree of rotational broadening. 
For the narrower lined stars (with WoSini ^ 150kms~^), a 
selection of metal lines and non-diffuse helium lines (Set 1 
in Table ^ was used. These had the advantage of being 
intrinsically narrow, whilst the former were not normally 
affected by nebula emission. For large rotational veloci- 
ties, reliable results could only be obtained from stronger 
absorption lines and in these cases a selection of diffuse 
and non-diffuse helium lines (Set 2 in Table [2]) was con- 
sidered. Some of the features are in fact close doublets or 
triplets, whilst additionally the helium diffuse lines have 
weak forbidden components (we note that blending of the 
He I 4026A feature with a Heii line should not be signifi- 
cant for our sample). For the former, the separations were 
less than the broadening due to rotation and indeed the 
instrumental profile, whilst the latter were only used for 
the more rapidly rotating stars where the profiles appear 
symmetric within the observational uncertainties. 

Normally it was not possible to obtain estimates from all 
the spectral lines in the adopted set for three main reasons. 
Firstly, some lines were too weak to provide reliable mea- 
surements. For example, the C ii and Mg ii features could 
not be used for the hotter stars, whilst the O ii and Si ill 
features were weak in the cooler objects. Secondly, blend- 
ing could affect some lines and in particular the He i line 
at 4121 A (blended with On) and the Mgn line at 448lA 
(blended with Aim). However the most serious constraint 
was the presence of strong nebular emission in the spectrum 
of some stars, which could not be reliably removed for our 
multi-fibre observations. This emission could affect the C n 
feature but more seriously the He i spectrum. The following 
approach was taken for increasing amounts of contamina- 
tion. 

1. The emission was snipped out of the spectra (see Fig. 
[1]), which was then regridded and analysed. The qual- 
ity of the Fourier transform, the consistency between 
the estimates from the FT and PF approaches and the 
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quality of the fit of the rotationally broadened profile to 
the observations was then used to assess the reliability 
of the estimate. 

2. The lines with the greatest nebular contamination (nor- 
mally the diffuse triplet transitions) were excluded. As 
can be seen from Table 4, there are no systematic dif- 
ferences between the estimates obtained from different 
lines. Hence the use of a subset of lines is unlikely to 
introduce significant biases. 

3. When a reliable estimate could not be obtained, an at- 
tempt was made to assign the star to a velocity range 
covering 40 kms^i (0-40, 40-80, 80-120 kms"^ etc.). 
The observed profiles were then compared with rota- 
tionally broadened profiles appropriate to this range to 
validate this decision. Hence although no reliable pro- 
jected rotational velocity could be estimated, it would 
be possible to include such stars in any statistical study. 

There were a number of narrow lined stars whose line 
width appeared to be dominated by the instrumental pro- 
file which corresponded to a FWHM of approximately 40 
kms~^. As would be expected, their Fourier Transforms did 
not show convincing minima with individual measurements 
showing a significant degree of scatter. Other stars showed 
a greater degree of broadening and yielded mean projected 
rotational velocities in the range of 30 to 40 kms~^. Here 
we have taken a conservative approach and have estimated 
rotational velocities only for stars with individual measure- 
ments from the Fourier Transform methodology that were 
consistently greater than 40 kms"-'^ and where the profile 
fitting implied a significant rotational broadening compo- 
nent. 

For each star, a projected rotational velocity was esti- 
mated only if there was at least three independent measure- 
ments. For the narrow lined stars, the average number of 
measurements was 6.9, with approximately 90% having five 
or more measurements, with the corresponding values for 
the broader lined stars being 6.3 and 93%. For the latter, an 
additional 16 stars were assigned projected rotational ve- 
locity ranges whilst one star, VFTS729, appeared to have 
asymmetric line profiles, suggesting a binary or composite 
nature, and was excluded. 

Tables 3 and 4 (only available online) lists all the in- 
dividual measurements using both the FT and the PF 
methodologies, together with the percentage depth of the 
feature as estimated from the profile fitting. The former is 
for the metal line set and the latter for the helium line set. 
In both tables stars that are of an uncertain status i.e. could 
possibly be binaries (see Sect. [5]) are identified with an as- 
terisk. The tables also give means and standard deviations 
with the FT estimates being weighted by the central depths. 
Using unweighted means would have lead to no systematic 
shift in the estimates with typical changes of 2% for line 
set 1, and a 1% decrease for line set 2 with typical changes 
again being 2%. Hence the choice of weighting is unlikely 
to be a significant source of error. We also list the mean 
values from the profile fitting; in this case the individual 
estimates were not weighted as these are better considered 
as providing upper bounds. Stars with narrow lined spectra 
that were instrumentally dominated have been assigned a 
pr ojected rotationa l veloci ty of 40kms~^. 

iTownsend eTaP (|2004[) and iFremat et all (|2005l) have 
investigated the effects of equatorial gravity darkening due 
to angular rotation (fi) on observed line profiles. They find 



Table 2. Sets of absorption lines used for projected rota- 
tional velocity estimates 



Set 1 Set 2 



Species 


A (A) 


Species 


A (A) 


He I 


4120 


He I 


4009 


He I 


4169 


He I 


4026 


Cii 


4267 


He I 


4120 


He I 


4437 


He I 


4143 


Mgii 


4481 


He I 


4387 


Sim 


4552 


He I 


4471 


Sim 


4567 


He I 


4713 


Sim 


4574 


He I 


4921 


On 


4661 






Hoi 


4713 







that for angular velocities near the critical limit (fie) the 
effects can be significant and lead to a systematic under- 
estima tion of the real projecte d rotational velocity. For ex- 
ample, 'Townsend et ajj (|2004[ ) considered a representative 
metal and helium diffuse line in stars with a spectral range 
B0-B3 (similar to that observed here) and for il/ilc=0.95 
and sinz — 1. They found the projected rotational ve- 
locity would be underestimated by between 9 and 22%. 
Unfortunately we do not know either the sin i values or at- 
mospheric parameters of our targets and are hence unable 

to systematically correct for this effect. 

Ho wever we can use the results of iMartavan et al.l 
(|2006D to obtain an insight into the magnitu de of any sys- 
tematic effects. They used the methodology ofl Fremat et al.l 
(2005) to correct observed projected rotational velocities 
for 47 Be-type stars towards the LMC cluster, NGC 2004. 
For ^l/ilc = 0.85 (the lowest value that they considered), 
the correction were 5.5% for stars with Vc sinz< 200kms~^, 
3.6% (200 <VcSini< 300kms-i), 3.3% (300 <VcSmi< 
400kms"i) and 3.0% (400 <VcSmi< 500kms"i). The sys- 
tematic errors in our estimates should be considerably 
smaller as it is unlikely that the majority of our targets 
are rotating with such high values of Q/Qc- For example, 
taking the median rotational velocity estimated from the 
de-convolution discuss ed in Sect. | 4 . 1 | and the critical veloc- 
ities from the grid of iBrott et al.l (|201lD for early B-type 
stars leads to a value of fi/J7c — 0.3 — 0.4. Hence we do 
not expect that there will be a significant effect on our 
sample as a whole, although a small number of individual 
estimates (where the target is rotating with near critical 
velocity) could be affected. 

3.2. Results and estimated uncertainties 

Before discussing the significance of the results presented 
in Tables 3 and 4, it is important to validate their re- 
liability. Several tests were undertaken as discussed be- 
low. Firstly one star (VFTS331) was inadvertently anal- 
ysed twice with excellent agreement; the mean of the FT 
estimates agreed to within 2kms~^ and that of the PF es- 
timates to within 3 kms^^. Secondly four targets in the 
overlap region between the narrow and broad lines stars 
(100 <,Ve s\ni< 200kms-i) were measured using both sets 
of lines and in Table [S] the mean values and the estimated 
errors are summarized. The estimates are in reasonable 
agreement with no evidence of systematic differences. The 
t-statistic and degrees of freedom for unequal sample sizes 
and unequal variance (Welch's t-test) are also listed. None 
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Table 5. Comparison of mean projected rotational veloci- 
ties (kms~^) deduced using the Fourier Transform method- 
ology from the two line sets. Also listed are the standard 
deviations and number (in brackets) of the estimates. The 
modulus of the t-statistic and number of degrees of freedom 
are also listed. 



Star 


Line set 1 


Line set 2 


\t\ 


df 


046 


169±15 (7) 


171±15 (7) 


0.2 


6 


060 


129±16 (9) 


129±8 (7) 


0.0 


6 


234 


157±12 (7) 


147±6 (7) 


L9 


6 


452 


125±15 (7) 


134±5 (8) 


L5 


7 



Narrow lined stars 




40 60 80 100 120 140 160 180 
FT estimates (km s"') 

Broad lined stars 




200 300 400 500 
FT estimates (km s"') 

Fig. 2. Comparison of projected rotational velocities esti- 
mated using the FT and PF methodologies. The upper fig- 
ure is for the narrower lined stars (line set 1), whilst the 
lower figure is for the broader lined (line set 2) sample. 
Also shown are the results of the simulation - crosses: He i 
4713A; asterisks: Sim 4552A or Hei 4387A. 



of these are significant at the 10% levels, but those for 
VFTS234 and VFTS452 are significant at at the 20% level. 

A comparison has also been made of the mean projected 
rotational velocities obtained from the FT and PF method- 
ologies. These are illustrated in Fig. [3] for the the narrow 
lined (line set 1) and broader lined (line set 2) subsamples. 
In general the results are in excellent agreement for both 
line sets. At lower projected rotational velocities (e.g. ^ 80 
kms~^), the profile fitting methodology yields systemati- 



cally higher estimates, which would be expected given that 
only rotational broadening was considered. 

To test this h ypothesis, we took the spectrum f rom a 
TLUSTY model flHubenvlll988t iHubenv et al.lll998l) . with 
an effective temperature of 25 000 K, a logarithm ic grav- 
ity o f 4.0 dex and an LM C composition (see Rvan s et al.l 
120031: iDufton et aDl2005L for details of the grid); This was 
convolved with a Gaussian profile to represent instrumental 
broadening and with different amounts of projected rota- 
tional broadening. The profiles were then analysed using 
the PF method to estimate the projected rotational veloc- 
ity. Three representative absorption features were consid- 
ered, viz. the Sim line at 4552A, the non-diffuse helium 
line at 4713A and the diffuse helium line at 4387A. These 
results are also shown in Fig. [2] and reproduce the trend 
in the differences between the observational estimates dis- 
cussed above. Indeed although they only represent a subset 
of lines for a single model, they provide strong evidence for 
the reliability of our Fourier Transform results. 

The sample standard deviations are typically 10% of 
the mean projected rotational velocities for both line sets. 
However the error in the mean values will be smaller and 
if the estimates are normally distributed, this would be the 
sample standard deviation divided by the square root of the 
number of observations. This would then lead to a typical 
uncertainty in the mean values of 4%. 

An independent error estimate can be found at higher 
projected rotational velocities, as here the profile fitting 
methodology should provide reliable estimates. We have 
therefore compared the two methodologies in these regimes. 
For line set 1, we considered all stars that had estimated 
projected rotational velocities of greater than 80 kms^^. 
For each star the ratio of the PF to the FT estimate was 
found, with the mean of these ratios being 1.01 ± 0.06. For 
line set 2 and limiting our sample to those stars with es- 
timates greater than 300 kms""'^ leads to a mean ratio of 
1.02 ± 0.04. Hence in these high projected rotational ve- 
locity regimes, the methodologies give excellent agreement, 
confirming the trends seen in Fig. [5] Additionally as the 
standard deviations for the ratios will include a contribu- 
tion from uncertainties in the PF approach, they are con- 
sistent with our estimated error of 4% for the means of the 
FT estimates. 

Rotational velocities have also been estimated for the 
apparently single 0-type stars in the Tarantula Survey 
(Ramirez- Agudelo et al., in preparation). They also used 
the FT technique but a different sets of lines reflecting 
the different temperature range. Due to uncertainties in 
the preliminary spectral classification, late O- and early B- 
type targets were included in both datasets (see Sect. [5]). 
Excluding targets that have only been assigned here to a 
velocity bin, there are 42 stars with independent estimates. 
There is no significant systematic difference between the 
two sets of measurements, with the mean value of the ab- 
solute differences being 3±15 kms~^and of the relative dif- 
ferences being 0±8%. For thirty seven of the targets, the 
difference is less than 10%, with the remaining five ranging 
from 10-30%. We discuss these discrepancy cases below: 

VFTS 226, 627: These targets have narrow hned spec- 
tra and are hence close to the resolution limit for getting 
reliable projected rotational velocities. Additionally macro- 
turbulence may be significant given their relatively large 
effective temperatures. Hence we do not consider the rela- 
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tively small absolute discrepancies (7-llkms ^) in the es- 
timates to be significant. 

VFTS 373, 543: Strong emission is present in both 
these spectra, making reliable estimates difBcult. Indeed 
our sample standard deviations are approximately 30 
km s"""^ (see Table 4), implying that the estimates are con- 
sistent virithin the uncertainties. 

VFTS 141: The spectrum of this star is well observed, 
with the nebular emission being weak. Although the two 
sets of measurements differ by 19%, when only lines that 
are common to both analyses are included this decreases to 
7%. Hence we believe that the measurements are consistent. 

3.3. Distribution of projected radial velocities 

A histogram of our projected rotational velocities, exclud- 
ing the targets with an uncertain status, is shown in Fig. [31 
A binsize of 40 kms~^ has been adopted as this is consis- 
tent with both our spectrally unresolved targets and those 
where only a velocity range was assigned. The distribution 
appears to be double peaked with approximately 25% of 
the sample, having i;oSini< 80kms~^. The sample should 
be relatively unbiased as it was selected on the basis of 
magnitude (with no colour cut-off) and then on the abil- 
ity to position the spectrograph fibres. The former could 
bias the sample towards targets with low extinction, whilst 
the latter could bias the sample away from regions with 
high star densities but neither of these would appear to 
be sufficient to explain the double peak. Another explana- 
tion of the bi-modal structure could be that many of the 
narrow-lined 'single' stars are in fact part of binary systems 
with synchronised orbits. We would expect such systems to 
have relatively small orbital periods with large radial ve- 
locity variations, whilst their radial velocity estimates will 
h ave a high accu r acy. A s such the radial velocity analysis 
of Duns tall et al.l (j2012l ) should have identified such targets 
and hence this bias should not exist. 

As discussed in Sect.^l we believe that our sample con- 
tains effectively all the targets with spectral types 09.7 or 
later but only partially samples those with spectral type 
09.5. Hence in Fig. [31 we also show a histogram with the 
09.5 stars excluded. We note that the two histograms are 
very similar with any differences being less than the random 
Poisson noise expected from finite sample sizes. Hence in 
the subsequent discussion of these data, we have considered 
the complete sample but excluded stars with an uncertain 
binary status. Our principle conclusions would remain un- 
changed if we had excluded the targets with a 09.5 spectral 
type and/or included those of uncertain status. 

3.4. Be-type stars 

The classical Be-type star is a B-type star that shows or 
has previously shown prominent emission features in its 
Balmer line spectrum, indicating the presenc e of a geo- 
mctrically fiattened, circumstellar disc (Ouirre nbach et al.l 
[1994L Il997t) . Although the mechanism for the Be phenom- 
ena remains unclear, such st ars are be lieved to rotate with 
velocities greater than 40% ()Cranme up to approxi- 

mately 90% of the critical limit, where the centri fugal force 
balances that of gravity (|Townsend et al.l [20041) . Reliable 
identification of Be-type stars in this sample was diffi- 
cult because of the strong nebular emission. Additionally 



sky subtraction did not work well for these fibre-fed spec- 
tra, because of the spatially highly variable nature of this 
emission. Hence to search for evidence of Be-type emis- 
sion, we examined red- region (HR15N) spectra processed 
without sky subtraction. Although there was no prospect 
of identifying weak or narrow stellar emission, broad Be- 
type Ha emission wings were identifiable in 52 out of 
the 289 constant-v elocity stars (including VFTS 102, see 
iDufton et al.ll20lj] ). and 8 out of 45 stars of uncertain status 
- i.e., about 18% of the sample. This fraction is comparable 
with canonical figures for LMC dust ers (e.g.. lMaeder et all 
119991: IWisniewski fc Biorkmanll2006[ ). even though we are 
only able to identify relatively strong emission. The Be-type 
candidates are identified in Tables 3 and 4 by crosses (+). 

In the bottom panel of Fig. [31 the projected rotational 
velocity distribution is shown for the sample excluding the 
Be-type candidates. As might be expected, this leads to a 
decrease in the number of rapidly rotating stars, whilst the 
bimodality is still present. There remains a small number of 
rapidly rotating stars (i;oSinz> 400 kms~^), although we 
cannot exclude them having weak Be-type emission. 

The median value of Vc sin i for the Be-star sample (ex- 
cluding VFTS 102) is 254 kms'^ (238 knis^i excluding the 
stars of uncertain status), and the mean 249 ±11 kms~^ 
(244 ± 12 kms~^). Although any bias introduced by the 
undetectability of weak stellar emission would presum- 
ably have skewed the results in favour of faster rota- 
tors, these values are somewhat smaller than the mean of 
^2 70 kms~^ for LMC fie ld and cluster Be stars reported 
bv lMartavan et al.l (|2006D . 

3.5. Comparison with previous studies 

Previous observational investigations of rotation in B-type 
stars have been discussed in the introduction. Here we limit 
the discussion to surveys in the Large Magellanic Cloud 
that are directly comparable with the current results. The 
FLAMES-I survey observed targets towards two LMC clus- 
ters. Nil and NGC2004 and provided estima tes of pro- 
jected rotational velocities ([Hunter et al.ll200R3) . These es- 
timates were not corrected for the systematic effect at near 
critical rotation discussed in Sect. 13.11 The rotatio nal ve- 
locity distribution inferred by [Hunter et al.l ([2008b[) had a 
modal value of 100km s~^corresponding to fl/flc — 0.2 and 
hence we would expect that most of their sample will not 
contain significant systematic effects. In Figs. [3| and [5l the 
projected rotational velocity distribution and cumulative 
distribution function are shown for the 92 targets that had 
no evidence for binarity, were not supergiants and had a 
spectral type of 09.5-B3 (consistent with our current sam- 
Ple |. 

[Martavan et al.[ ([20061 ) provided estimates for 106 B- 
type and 47 Be-type stars towards NGC2004. Supergiants 
and spectroscopic binaries were also excluded in this study, 
and most investigated stars had masses in the range 4 M0to 
20 Mq. Again limiting the sample to the 146 09.5-B3 leads 
to the distributions shown in Fig. [H and [51 

Neither of these tw o studies iden t ified a bi-modal dis- 
tribution. For example, [Hunter et al.[ ([2008btl fitted a sim- 
ilar histogram (but with a different bin size and sample; 
their Fig. 9a) by a uni-modal distribution. However tailor- 
ing their sample to match our range of spectral types leads 
to evidence for a bi-modal distribution, although this must 
be considered as marginal given their small sample size. By 
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Fig. 3. Histograms of the estimated projected rotational 
velocities, excluding stars of uncertain status. Top: all tar- 
gets; Middle: excluding 09.5-type targets; Bottom: exclud- 
ing Be-type candidates 



contrast, the larger sample of iMartavan et al] (|2006D shows 
no clear evidence for bi-modality, although this cannot be 
excluded. 

To investigate whether these samples originated from 
the same parent population as our current sample, we 
have used the two sample Kolmogorov-Smirn qv and Kuiper 
tests. T hese implied that for the sample of iHunter et al.l 
(|2008b[) . this was unlikely (P r^ 1 x 10"^ and 6 x 10"^ 
respectively) but for that of IMartavan et ahl (|2006[ ). the 
probability was relatively high (P ~ 0.4 and 0.2 ). This 
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Fig. 4. Histograms of the esti mated pro.iec t ed rota tional 
velocities for the sur veys of IHunter et all ()2008b[ ) and 
IMartavan eraD (|2006D . 



is consistent with cumulative distribution functions plotted 
in Fig. [5j 

These comparisons must be treated with caution as the 
samples remain poorly matched. Firstly the selection cri- 
teria adopted for the three surveys will have introduced 
differen t observationa l biases . For example both the sam- 
ples of IHunter et aP (|2008bD and IMartavan etal! (|2006[) 
contained evolved main sequence stars (see their Figs. 7 
and 5, respectively). As the 30 Doradus region is thought 
to be young, our sample should contain a higher fraction of 
relatively unevolved B-type stars, although we do not cur- 
rently have spectral types or luminosity estimates to con- 
firm this. Secondly the previous studies observed a mix- 
ture of cluster an d field stars, wi t h the latter dominat- 
ing. For example, IMartavan et al.l (|2006l ) estimated that 
r-^75% of their sample were field stars; given the complex 
nature of the 30 Dor region, the fraction of field stars (in 
the traditional sense) in the curre nt sample will be rela - 
tively low. Thirdly the sample of IMartavan et al.l (|2006|l 
was obtained at two epochs and they only identified ap- 
proximately 10% of their sample as b i naries. By contrast 
for the sample of lHunter et"all (|2008bD . lEvans et al.l (|2006D 
had identified 23% of the NGC2004 targets and 36% of 
the Nil targets as binaries and these have been excluded 
from the sample considered here. Hence given that the 
Tarantul a survey has a lower limit of 37% on the binary 
fraction (jPunstall et al.l[2012[ ). we conclude that the sam- 
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Fig. 5. Cumulative distribution function for the our esti- 
mated projected rotationa l velocities (so li d) , tog ether with 
thos e for the surveys of Hunter et al.l ()2008b() (dashed) 
and iMartavan et al.l ()2006[) (dot-dashed) . For the current 
dataset, the unresolved spectra have been given an arbi- 
trary UeSini of 20 kms~^ 

pie of lMartavan et all (120061 ) probably contains a significant 
fraction of binaries. 

4. Rotational velocity distribution 

The relatively large size of our sample has been exploited to 
investigate its intrinsic rotational velocity distribution. We 
first describe our methodology, and then outline our results 
and investigate their robustness. 

4.1. Deconvolution of the Vc sin i distribution 

If the stellar rotation axes are randomly distributed, it 
is possible to infer the probability density for the rota- 
tional velocity distribution (P(i;o)) from that for WoSini. 
Previously several authors have investigated this by assum- 
ing an analytical function for P(we) and then optimising the 
fit with observation after convolution to allow for random 
inclinations. For example, for the analysis of the FLAMES- 
I datasets, rect angular, Max wcllian (Mokicni ct al. 2006|) 
and Gaussian (iDufton et al.l [2006bl: iMokiem et al.l \26m 
iHunter et al.ll2008b[ ) distributions were adopted. Given the 
structure in the projected rotational velocity distribution 
discussed above, it is not clear what adopted function would 
be appropriate. 

We have therefore attempted to deconvolve the pro- 
jected rotational velocity distribution to directly infer that 
of the rotatio nal ve l ocitie s. We have adopted the iterative 
procedure of 'Lucy' fl974) which if it converges will con- 
verge to the point of maximum likelihood. For the continu- 
ous function, representing the observed projected rotational 
velocities {4> in the notation of Lucy), we followed his recom- 
mended procedure. The initial estimate for the rotational 
velocity distribution (ip'^ in the same notation) was taken 
to be a Gaussian with a positive abscissa; tests showed that 
this choice did not affect t he fina l iterat ive solution, consis- 
tent with the discussion of iLucvl (|1974[ ). 

The effectiveness of this procedure has been tested 
by taking simulated distributions for P(tie) and then us- 
ing a Monte- Carlo technique to generate P(i7eSini). These 



were then de-convoluted to obtain an estimate for P(fc). 
Simulated distributions considered included Gaussian, tri- 
angular and step functions. In general the deconvolution 
retrieved the large scale structure, while the estimated pro- 
jected rotational velocity distribution, ((f> in Lucy's nota- 
tion), agreed well with the simulated P(fo sini). Additional 
small scale structure was present especially where the 
adopted rotational velocity distribution was discontinuous. 
Also structure varied from simulation to simulation and ap- 
peared to be connected with that due to the Poisson noise. 
The only significant concern to arise in these simulations 
was that on occasions a small excess of slowly rotating stars 
were found even though there was excellent agreement be- 
tween the simulated projected rotational velocity distribu- 
tion, P(woSinz), and that inferred from the deconvolution 

procedure (^). In turn this implied that relatively small 
errors in the observed projected rotational velocity distri- 
bution could propagate through to that for the rotational 
velocities, P(i;c). 

Our adopted pro cedure is similar to that used by 
IZorec fc Roved (j2012L and references therein) to investi- 
gate the rotational velocity distrib utions for Gala c tic A- 
type stars. However we differ from IZorec fc Roved ((20TI 
in not preprocessing our observational projected rotational 
distributions using a Gaussian filter. Our approach will not 
remove spurious structure in our inferred rotational veloc- 
ity distribution due to the Poisson noise in the observa- 
tional distributions. However it will also illustrate the lim- 
itations of this approach and its sensitivity to the assump- 
tions made. 

The choices of the number of iterations and the bin 
size for the observed projected rotational velocity are im- 
portant. For the former, the iteration initially converges 
rapidly and then tends to accentuat e small scale structure 
that may not be real. iLucvl ()1974[ ) in his two rotational 
velocity examples only considered two iterations. In Fig. 
El we show (for an observation bin size of 40 kms"""^), the 
originally assumed distribution of the rotational velocities 
and the solutions after two, four, six and eight iterations. 
The solution after two iterations shows the major features 
of the distribution with further iterations introducing addi- 
tional structure at approximately at 70, 270 and 420kms~^. 
Although the signal will not be band-limited, this struc- 
ture would appear to contain frequencies higher than the 
Nyquist frequency (or equivalently in the rotational veloc- 
ity domain smaller than half the bin size) . Additionally the 
errors due to Poisson statistics map onto a typical uncer- 
tainty of approximately 5xl0~^ in P{vc) and hence this 
additional structure is probably not significant. 

In the case of the adopted observational bin size it is 
necessary to balance the resolution in the velocity domain 
against the (Poisson) noise in the star counts. We have 
taken a pragmatic approach by considering two different 
possible bin sizes of 20 and 40 kms~^, with for the former 
the spectrally unresolved targets being equally distributed 
between the first two bins. The results are illustrated in 
Figs. [6] and [7] for different number of iterations. The large 
scale structure of the estimates of the rotational velocity 
attribution are similar but the smaller bin size leads to more 
small scale structure. Although this will to some extent 
refiect the higher sampling rate, we believe that mainly 
that it arises from random errors in the Ve sini statistics. 
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Fig. 6. Distribution of de-projected rotational velocities 
from our iterative procedure for the observed Ve sin i distri- 
bution binned to 40 kms^^. The initially assumed Gaussian 
distribution (solid line) and the results after two (dashed), 
four (dot-dash), six (double dot-dashed) and eight (solid) 
iterations are shown for Vc< 500 kms. 
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Fig. 7. Distribution of de-projected rotational velocities 
from our iterative procedure for the observed Vc sin i dis- 
tribution binned to 20 kms^"'^. Notation is identical to Fig. 

m 



On the basis of the above tests, we have adopted the 
larger bin size together with 4 iterations in our analysis and 
discussion given below. We note that the major conclusions 
would not have been significantly affected by other choices. 

4.2. Results 

In order to illustrate our results. Fig. [5] shows the observed 
Ve sini histogram from Fig. [31 normalised to unity when in- 
tegrated over velocity. Also shown are the estimates after 
four iterations of the probability density of the projected 
rotational velocity distribution (0* in Lucy's notation) and 
of the corresponding probability density for the rotational 
velocity, P(uo) ('0'* in Lucy's notation). There is good agree- 
ment between the observed Vc sin i distribution and that 
inferred (0'*) from the iterative process. Additionally the 
estimate of the rotational velocity distribution appears rea- 
sonable with the large scale structure of the sin i distri- 
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Fig. 8. Histogram of observed normalised Vc sini distribu- 
tion binned to 40 km s~^ . Also shown are the estimates after 
4 iterations of the probability densities for the projected ro- 
tational velocity distribution, P(wcSini) (0^ in Lucy's no- 
tation - dot-dashed line) and for the rotational velocity, 
P(uc) in Lucy's notation- solid line). 



butions being conserved but moved to larger velocities due 
to the effect of inclination. 

In Table [5] our best estimates of the probability distri- 
bution function are tabulated. For smaller rotational veloc- 
ities, the estimates have been rounded to SxlO^'"^ but for 
larger velocities a small interval was chosen. We emphasis 
that these values are subjective as they require a judge- 
ment about the validity of the smaller scale structure. Also 
listed in this table is the cumulative distribution function 
(cdf) corresponding to this best estimate and because this 
is truncated at Vc< 600 km s^^, it does not reach unity. We 
recommend that: 

1. either the range up to the critical velocity is populated 

2. or that the cdf and estimate of P(wc) are renormalised 

depending on the nature of the application. 

All the de-convolutions show a double peak in the 
equatorial velocity distribution consistent with that ob- 
served in the projected rotational velocity distribution. 
Approximately one quarter of the sample have a rotational 
velocity of less than 100 kms~^ and there appears to be 
another peak 250-300 kms~^ although this is complicated 
by the varying degree of small scale structure found in the 
different de-convolutions. 

We have attempted to fit the de-convolved distribution 
using analytical functions. Initially two Gaussian profiles 
were adopted but these were found to give a relatively poor 
fit to our estimated distribution. Two M axwellian func- 
tions as discussed by IZorec fc Roved (|20T1 were also con- 
sidered but would not have reproduced the significant value 
of P{vc) as 0. We therefore recommend that either the 
values listed in Table [5] or analytical fits appropriate to the 
specific applications are adopted. 

The most distinctive feature of the Vc sin i distribution 
is its bi-modal nature, which is also present in the estima- 
tion of the rotational velocity distribution, illustrated in 
Fig. |8l Although our sample contains approximately 300 
targets, it will still be subject to significant random sam- 
pling errors due to its finite size. In Fig. |21 the population 
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Table 6. Estimates of the probability density of the rota- 
tional velocity, P(wc), and its cumulative distribution func- 
tion, cdf. 



Vc 


P(t'c)xlO^ 


cdf 





2.20 





20 


2.40 


0.046 


40 


2.45 


0.094 


60 


2.55 


0.144 


80 


2.25 


0.192 


100 


1.40 


0.239 


120 


1.00 


0.253 


140 


0.75 


0.270 


160 


1.05 


0.288 


180 


2.40 


0.322 


200 


3.00 


0.377 


220 


2.85 


0.435 


240 


2.75 


0.492 


260 


2.90 


0.548 


280 


3.20 


0.609 


300 


3.30 


0.674 


320 


3.20 


0.739 


340 


2.50 


0.796 


360 


2.00 


0.841 


380 


1.80 


0.879 


400 


1.50 


0.912 


420 


1.10 


0.938 


440 


0.75 


0.956 


460 


0.70 


0.971 


480 


0.50 


0.983 


500 


0.21 


0.990 


520 


0.09 


0.993 


540 


0.08 


0.995 


560 


0.06 


0.996 


580 


0.04 


0.997 



of the two 40 kms~^ bins with the lowest projected ro- 
tational velocities have values of approximately 40 corre- 
sponding to an estimated standard error of approximately 
6. Hence even decreasing the populations of both these 
bins by twice this estimated error (note that if these bins 
had been overpopulated our estimate of the standard error 
would also be too high) would not remove this bi-modal be- 
haviour. Additionally Kolmogorov-Smirnov tests using uni- 
modal distributions (for example a single Gaussian fit to 
our estimated rotational velocity distribution) lead to very 
small probabilities that they were the parent populations. 
Hence we conclude that the projected rotational velocity 
distribution (and hence the underlying rotational velocity 
distribution) is bi-modal, although we accept that our sam- 
ple size limits the information on the structure of the two 
components. 

In Sect. 13.51 we discussed the projected rotational ve- 
locities in two other LMC samples. Those of iHunt er et al.l 
(|2008bD showed some evidence for a bimodal distribution 
but that of iMartavan et al.l (|2006[ ) appeared unimodal. 
However given the differences in the samples in terms of 
age, fraction of field stars and undetected binaries, we do 
not consider this discrepancy to be significant. 

5. Origin of bi-modal distribution 

The bi-modal distribution of rotational velocities implies 
that our stellar sample is composed of, at least, two dif- 



ferent components. Which physical processes lead to the 
existence of these two components? Why do most B-type 
main sequence stars rotate fast, while about a quarter or 
so rotate slowly? Whatever the answers to these questions, 
they are not contained in the standard evolutionary theo- 
ries of single stars, which predicts only small changes of the 
surface rotational velocity during core hydrogen burning 
(see, for example. iBrott et al.ll201ll ). 

The two components, which we find, could either corre- 
spond to different ages, different star formation conditions, 
or emerge as a consequence of differences in the evolution of 
stars (e.g., in close binaries). We emphasize again that age 
differences can only lead to the difference in the mean rota- 
tion rate of the two components unless the standard stellar 
evolution picture is wrong or incomplete. Differences in age 
and star formation conditions could produce stellar com- 
ponents with different spatial or kinematic properties. We 
inspect our data in this respect in the next two subsections. 
Thereafter, we discuss the possible impact of non-standard 
stellar evolution. 

5.1. Spatial variations 

The 30 Doradus region has had a complex star-formation 
history over the past few tens of Myrs, as demon- 
strated by the fiv e disti nct stellar groups identified by 
IWalborn fc Blades! ()19970 . For instance, the VFTS ob- 
servations include targets in the cluster Hodge 301, 
which is ^ 20Myr older than the central parts of 
30 Doradus (|Grebel fc Chul |2000[ ). Indeed, even in the 
centra l 30 Doradus cluster (NGC2070), ISelman et all 
(1 1 9991) identified three distinct temporal components, 
with Sabbi et al.l ()2012|) recently identifying two spatially- 
distinct components which appear to have an age differ- 
ence of '^IMyr. An obvious test is thus to check whether 
the slow rotators concentrate in particular regions of the 
VFTS field. In Fig [HI the spatial positions of the sample 
of our targets are illustrated, together with their estimated 
projected rotational velocity. The concentrations of stars 
at specific positions follows directly from the morphology 
of the region. 

From this figure, there appears to be no evidence for any 
spatial segregation with projected rotational velocity. For 
example, the percentage of targets with values of Vc sin i less 
than lOOkms"^ within 0.1 degrees of the central concentra- 
tion of stars associated with Hodge 301 is 33%, in excellent 
agreement with the 33% found for the whole sample. 

To further assess any spatial variation, the sample was 
split into five velocity bins, each of 100km s~^, except 
the top bin which encompassed all objects with UoSini> 
400km s~^. Cumulative radial-distribution functions were 
calculated for each bin. Two-sample Kolmogorov-Smirnov 
(KS) tests were then performed on all possible pairs of 
the radial-distribution functions. No statistically significant 
variation was found between the pairings, suggesting the 
spatial distribution does not influence the sin i distribu- 
tion. 

I n this context it i s interesting to consider the results 
from lWolff et al.l (|2008l) for projected rotational velocities of 
24 OB-type stars in R136. Their sample stars were located 
within ~25" of the centre of R136 (equivalent to a projected 
distance of ~6pc), i.e. the majority of their stars lie within 
the dense central part of the FLAMES field that was not 
observed with the Medusa fibres. From comparisons with 
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Fig. 9. Spatial position of our targets with different colours 
used to identify different ranges of projected rotational ve- 
locity, according to the colour bar at the right side. The 
relatively old cluster Hodge 301 lies at a right ascension of 
84.50° and a declination of -69.06°. 



Table 7. Average radial velocities for sub-samples with 
vsini smaller and larger than lOOkm/s (cf. Fig [3]), com- 
pared to the full sample. 



Sample 


Number 


(iJrad) kmS ^ 


feSinj< lOOkms"^ 


89 


269±17 


VcSmi> 100 kms^^ 


172 


267±18 


All 


261 


267±17 



In Table [71 we summarize the mean and standard de- 
viations for the different samples, together with the sam- 
ple sizes. Although reliable values were not available for 
approximately 10% of the targets, we do not expect this 
to introduce a significant bias. For the two samples, both 
the means and standard deviations are effectively identical. 
Coupled with the results presented in Sect. 15. ll we conclude 
that there is no evidence that the bi-modality arises from 
different populations, due to either sequential star forma- 
tion or from contamination with field stars. 



results for other samples in the LMC (including those from 
Hunter et al. 2008b), Wolff et al. concluded that the mean 
rotational velocity of the massive stars in R136 is larger 
than that found from other samples in the LMC, both for 
stars in (less dense) clusters and in the field population. 

IWolff et al.l (|2008D noted that R136 lacks the popula- 
tion of very slowly-rotating stars seen elsewhere (and which 
also features in the results presented here). They divided 
their sample into high- and low-mass bins, with inferred 
evolutionary masses in the range of 15 < < 30 Mq and 
6 < Af, < 12 Mq, respectively. The majority of the stars in 
the lower-mass bin lie below the faint magnitude cut-off for 
the VFTS, but the lack of slow rotators is also evident in 
their high- mass sub-sample (see their Fig. 11) - they ar- 
gued that these results might refiect the initial conditions 
within which the stars formed, i.e. the relationship of rota- 
tional velocity with stellar density. As noted above, we see 
no compelling spatial trends in the results from our signifi- 
cantly larger VFTS sample, with the caveat that our results 
do not include stars in the central few pc of R136. 

5.2. Radial velocity variations 

A bi-m odality found in Galactic targets by iHuang et al.l 
(|2010t) has been interpreted as due to a field star and a 
cluster population. Although we see no evidence for this in 
the spatial distribution of targets, we have also investigated 
the estimated stellar radial velocities of our sample. These 
have been measured Kennedy et al. (in prep.) using pro- 
file fitting to selected spectra l features w i th a m ethodology 
similar to that discussed by ISana et alj ([20123) . We have 
divided our sample into two groups with projected rota- 
tional velocities less than or greater than 100 kms~^. The 
former will be a mixture of slowly rotating stars and those 
viewed at low angles of inclination; however the contami- 
nation should be relatively small given that the mean value 
of sini for random inclinations is approximately 0.78. The 
latter will be stars with rotational velocities in excess of 
100 kms-^ 



5.3. Binary evolution 

The fraction of G alactic massive clo se binaries is high, 
i.e., at least 50% (|Sana et al.ll2012bD . Similar high frac- 
tions have been found f or the O- and B-type VFTS sam- 
ples ISana et all (|2012aL cf.. Sect. [5]). Since the time evo- 
lution of the rotation of massive main sequence stars can 
be strongly a ffected by the interaction with a companion 
(jLangeil 2012fl . this could in principle affect the distribu- 
tion of rotation rates of our sample. The reason is that, 
while radial velocity variab l es have been eliminated from 
our sample, Ide Mink etHI (|201ll l2012bD found that the 
vast majority of post-interaction objects are either single 
stars or appear as such. 

De Mink et al. (2012b) showed that for constant star 
formation, the fraction of apparently single post-binary in- 
teraction objects can be as high as 30%. This fraction would 
be considerably larger in a sample in which radial veloc- 
ity variables are removed, as except for a small fraction 
of Algol-type systems, those are almost exclusively pre- 
interaction binaries. For the current sample, for which the 
age distribution is not yet known, on one hand the frac- 
tion of apparently single close binary products is expected 
to be smaller than 30% due to the expected presence of 
early B-type stars which are coeval with the core region of 
30 Doradus. On the other hand, this fraction is enhanced 
due to the removal of the radial velocity variables from our 
sample. In conclusion, at present a substantial contribution 
of stars that are the product of binary evolution, can not 
be excluded. 

The bina r y po pulation synthesis models of 
Ide Mink et al.l |20123) imply that binary interactions, 
i.e., mass and angular momentum transfer and stellar 
mergers, produce a population of fast rotators, with 
rotational velocities up to about 600 km/s. In principle, 
binary interaction can also produce slow rotators, for 
example , as a result of spin down by tides. However, 
dc Min k et al.l (|2012aD find that the effect of this on a 
stellar population is very small. The low velocity peak of 
our distribution (see Fig. 8) contains about 25%, which is 
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too many to be explained solely as the products of binary 
evolution. Indeed, to produce the low velocity component 
of our distribution by close binary effects could require 
the hypothesis that mass transfer and/or mergers lead to 
significant magnetic fields in the binary re mnant (cf., Sect. 
5.4) as supported by the recent discovery bv lGrunhut et al.l 
(2012b') of a strong magnetic field in the companion of 
Plaskett's star. 



5.4. Magnetic fields 



IZorec fc Roveij (|2012l . and references therein) have under- 
taken an extensive study of the rotational velocities in 
Galactic late B-type and A-type stars, using a similar 
methodology to estimate rotational velocity distributions. 
Their large dataset allowed them to study variations as a 
function of age and mass. They find that stars with masses 
^ 2.5 M0 have a uni-modal distribution of rotational ve- 
locities, whilst the more massive stars (2 .5-3.9 Mfr^) show a 
bi- modal structure. IZorec &: Roveij ()2012[ ) fit the latter with 
two Maxwellian functions and find that the percentage of 
targets in the low velocity Maxwellian varies from 8-20% 
depending on the mass range considered. Additionally they 
find that the peak of the higher velocity Maxwellian lies in 
the range 200^t;o^240kms"^ for the highest masses they 
consider. 

These results ar e not directly compar able with those 
presented here since IZorec &: Roveij (|201^ sampled differ- 
ent field and cluster environments. Additionally our sam- 
ple is not large enough to be subdivided into smaller mass 
ranges, while we do not have age estimates for our targets. 
However, both the fraction of slow rotators (25% for our fit; 
8 — 20% in the higher mass samples of [Zorec fc Rover ) and 
the peak value of the high velocity component (280 kms^^ 
an d 200-240 km s~^ resp ectively) are similar. 

I Zorec fc Roveij (|2012D discuss two possible explanations 
for the bi- modality in their distributions. The first is that a 
fraction of stars loses angular momentum via tidal breaking 
within a bin ary system. For our sam ple, we can exclude this 
scenario, as Ide Mink et al.l ()2012bf ) show that the number 
of affected st ars is much tog small. The second possibility 
discussed by IZorec &: Roved (|2012f ) is that all stars were 
born rotating rapidly, but that a fraction of them had or 
acquired magnetic fields. Indeed, the fraction of magnetic 
main sequence sta rs in the mass range 2.5-3. 9M f:iis believed 
to be about 15% (jPonati &: Landstreetll2009l ). whilst their 
rotation may be slowed down by magnetic breaking. 

It has become evident in recent years that magnetic 
fields could potentially play a significant role in the evo- 
lution of some massive stars. In the LMC sample of earl 
B-type stars discussed in Sect. 13.51 iHunter et al.l (|2008. 
found a sub-sample of 15% of their stars with slow ro- 
tation and high nitrogen enrichment. While direct mag- 
netic field m e asure ments in these stars are no t yet feasible, 
iMorel et al.l (|2008D and iHubrig et ahl (|201ll) have identi- 
fied Galactic analogues for which the incidence of magnetic 
fields is indeed high. Grunhut et al. (2012a) find magnetic 
fiel ds in ^10% of their e arly B stars. 

iMevnet et all (|201lD and iPotter et all (|2012ajlb[ ) have 
explored stellar evolution models for magnetic massive 
m ain sequence stars. Based on a model for magnetic braking 
of lud-Doula fc Owockij(|2002j ). they both find that magnetic 
massive main sequence models can in fact s pin- down for 
reasonable field strength. iPotter et al.l ([2012^ performed a 



p opulation synthesis m atched to the B-type stellar sample 
of lHunter et al.l ([2008a) . They produced a bi- modal velocity 
distribution, with a significant component of slow rotators 
(ve < 60km/s). However their distribution of fast rota- 
tors peaks at Ve — lOOkm/s (possibly due to their choice 
of the initial rotational velocity distribution) , compared to 
Ve — 250km/s for our sample (Fig. |S]). 

Of course, if magnetic spin down was responsible for 
the low-velocity component in our rotational velocity dis- 
tribution, one may wonder about the origin of the mag- 
netic fields in these main sequence stars. Three possi- 
bilities are being discussed in the literature. The fields 
could be fossil, i.e. inherited from the star formation pro- 
cess (jPonati fc Landstreetl I2OO9I) . Alternatively, the fields 
could be generated by strong binary interaction, i.e. bi- 
nary merger, mass transfer or co mmon evelope evolution 
(jFerrario et al.ll2009l:lLangeij|2012D . Finally, the fields could 
be generated by a dynamo process inside the the main 
sequence star (ISpruitJ I2OO2I: IPotter et al.ll2012bh . We em- 
phasize that these three alternatives predict different time- 
dependances of the stellar rotation and magnetic field in- 
cidence, which might be testable if the age distribution of 
our sample stars can be established. 

An additional important braking mechanism could be 
mass loss via stellar winds. Winds are thought to be ubiq- 
uitous amongst massive stars, and despite the spectroscopic 
evidence for high-speed outfiows being most prevalant 
amongst the most luminous 0-type objects (in the form 
of P Cygni ultraviolet lines), they are also observed in 
the lower mass B-type stars. Although mass-loss rates are 
higher in O stars, wind braking (Langer 1.998) could be 
relevant for B-type stars given their longer evolutionary 
timescale s. One interesting possibility is that of bi-stability 
braking (jVink et al.l l2010l) , which might lead to a popu- 
lation of both slow and fast rotators, and hence play a 
role in the observed bi-modality. This could be tested when 
the stellar parameters of our sample stars have been estab- 
lished. 



6. Conclusions 

From the VFTS B-type stellar sample (approximately 540 
stars), we have constructed a sample 289 early B-type stars 
(plus 45 stars of uncertain binary status) by removing ra- 
dial velocity variable stars and supergiants. These will have 
masses in the range 6 to 16 Moand are centered on the 30 
Doradus region in the LMC. We have estimated their pro- 
jected rotational velocities, and find that their distribution 
shows a distinct bi-modality, which neither correlates with 
position in the field nor with radial velocity. 

We have deconvolved the Vc sin i distribution to ob- 
tain the most likely underlying distribution of intrinsic ro- 
tational velocities. This shows a low velocity component 
containing about one quarter of the sample stars, with 
Vc <100km/s, and a second component of stars rotat- 
ing predominantly with velocities in the range 200km/s< 
Vc < 350km/s. 

This distribution is not inconsistent with previous LMC 
surveys, despite a bi-modality having not been clearly 
demonstrated before. A conclusive interpretation of our 
finding is not yet possible, in particular because the age 
distribution and chemical composition of the sample stars 
are still unknown. An origin of the bi-modality due to evo- 
lutionary effects relating to binarity and/or magnetic fields 
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or that the velocity distribution was inherited in the star 
formation process remain possibilities. 
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